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Chemical abundances in M31 from H II regions 
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■ ABSTRACT 

We have obtained multi-slit spectroscopic observations from 3700 A to 9200 A with 
i— i. LRIS at the Keck I telescope for 31 H Ii regions in the disk of the Andromeda galaxy 

■ (M31), spanning a range in galactocentric distance from 3.9 kpc to 16.1 kpc. In 9 Hn 
\ regions we measure one or several auroral lines ([Oin]A4363, [Nll]A5755, [Sm]A6312, 

[On]A7325), from which we determine the electron temperature (T e ) of the gas and 
derive chemical abundances using the direct T e -based method. We analyze, for the first 
time in M31, abundance trends with galactocentric radius from the direct method, and 
O 1 find that the Ne/O, Ar/O, N/O and S/O abundance ratios are consistent with a con- 

s'^ | stant value across the M31 disc, while the O/H abundance ratio shows a weak gradient. 

r/} . We have combined our data with all spectroscopic observations of Hn regions in M31 

available in the literature, yielding a sample of 85 H n regions spanning distances from 
3.9 kpc to 24.7 kpc (0.19-1.2 R25) from the galaxy center. We have tested a number 
. of empirical calibrations of strong emission line ratios. We find that the slope of the 

' oxygen abundance gradient in M31 is —0.023 ± 0.002 dex kpc -1 , and that the central 

I oxygen abundance is in the range 12+log(0/H) = 8.71 -8.91 dex (i.e. between 1.05 

and 1.66 times the solar value, for 12+log(0/H)©=8.69), depending on the calibra- 
l/"") ' tion adopted. The H 11 region oxygen abundances are compared with the results from 

other metallicity indicators (supergiant stars and planetary nebulae) . The comparison 

■ shows that Hn region O/H abundances are systematically ^0.3 dex below the stellar 
' ones. This discrepancy is discussed in terms of oxygen depletion onto dust grains and 

£Nl ■ possible biases affecting T e -based oxygen abundances at high metallicity. 
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1 INTRODUCTION iDavidge etall 120121 : iTanaka et ail l20ld) which have radi- 
cally changed our knowledge on the formation and evolution 

M31, the Andromeda galaxy, is our nearest neighbour and „ , , 01 T , . , , . , ,. , , „., 

J of M31. It is now recognized that interactions and satellite 

the largest galaxy of the Local Group and these properties , . , , . , , ., n c ., , , . c 

f 1 » , f r , accretions have been widely responsible tor the evolution 01 

make this galaxy a fundamental benchmark tor understand- , , , , , , , | T , , — — TNT-T^TTTI [7! — — T"^ — 71 

the Andromeda galaxy (e.g. ibata ct al. 2001; Bernard et al. 

mg spirals and their evolution. M31 s proximity allows tor |„„ ln | rf: — rr— — 7— nioniol It T 1 llnnr,J\ mL 

1 \ 2012; Davidgc et al. 2012: Hammer et al. 2007). The mter- 

spatially detailed studies of its constituent stellar popula- , ,, ,. , , . r ,. 

J f stellar medium also shows signatures ol disruptions possi- 

tions, even though the large angular extent on the sky ham- nj . ,, , . , , |Vt ■ id^ool 

' . . ™ , r , n i 1 bly resulting through these interactions (e.g. Unwml 119831 : 

pers a global view of this galaxy. The advent of large field | Co rbelli et aljB . The cu rrent star formation rate (SFR) 

cameras and multi-object spectrographs has been critical ,, » _i j,, r . nl . — n[7-T^7T^l\ -n ,1 

P P across the disc, ~1 M yr (Williams 2003), is significantly 

for a growing amount of information on important prop- , ,, . , , ,. 

lower than m the Milky Way and most of the star-tormmg 

erties of M31. These improvements m observational capa- ., ,. , . , . . . „ , , 

sites are distributed m a prominent ring ot approximately 

bilities have permitted, for example, the discovery of a gi- , , , . ,. . n , , , 

. \ \ '17; -11 " 1, , ~10 kpc galactocentric radius, this morphology has been 

ant stellar tidal stream m the halo Ibata ct al. 2001) and ,. . , ... , . , , , r ,, , . , 

., , „ , , „' : — . — ~ , . predicted with dynamical models ot the recent interaction 

detailed studies ot the resolved stellar populations both m ,,,„, .,, . , , , rr^j — ; — 7 — ri 

1— t — <- 1 1 -1 1 1 ot M31 with a companion galaxy possibly M32, Block ct al. 

the disc and in the halo (e.g. iFerguson et all 12003 . 12003 : r^b. v 



The chemical evolution of M31 has also been ex- 
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2006). 
T 

tensively stucuea m comparison witn tne MiiKy way 
l|Marcon-Uchida. Matteucci. fc Costall2010l : lYin et alj|200a : 



2 A. Zurita and F. Bresolin 



iRenda et al.ll2005h . Both Local Group galaxies seem to share 
common features in their chemical evolution, but in order to 
explain the observations in detail, models require different 
formation histories for the two galaxies. In particular, M31 
must have been more active an d efficient in the past in form- 
ing s tars than the Milky Way ()Yin et al.ll2009l ; IRenda et all 
l2005h . 

The present-day chemical element abundances and 
their spatial distribution across a galaxy are key observa- 
tional properties f or constraining different mo dels of chem- 
ical evolution (e.g. iMarcon- Uchida et ai]|201Ch . Hn regions 
have been traditionally the main targets used to trace 
the present-day elemental abundances across galaxy disks. 
However, nebular metallicity data on M3I are surprisingly 
sparse. Since the firs t H H region emission-line strength mea - 
surements made by iRubin. Krishna Kumar, fc Fordl (|l972l ^) 
in the early 1970's, in which the existence of an abun- 
dance gradient was suggested, only a handful of authors 
have tried to meas u re th e abundance gradient of M31. 
iDennefeld fc Kunthl (| 19811 s ! analyze a sample of 12 su- 
pernova remnants (SNR) and 8 Hn regions. They re- 
port gradients in the N/H and N/O ratios across the 
galaxy disk, but no gradient in the distribution of oxy- 
gen abundance with galactocentric radius, within the large 
uncertainties in t heir abundances determinat i ons. In a 
subsequent study, iBlair. Kirshner. fc Chevalier! l|l982f ) ob- 
tain an oxygen abundance gradient with a slope of ap- 
proximately —0.027 dex kpc _1 and a super-solar oxygen 
abundance in the M31 center (12+log(0/H) > 9.0) from 
the application of empirical strong-line nebular metallic- 
ity determination methods to the observed line fluxes in 
a sample of 11 Hn regions. In the same paper they 
calculate the oxygen abundance gradient from the spec- 
tra of 9 SNRs, and these show no significant gradient 
across the disk, with a mean abundance 12+log(0/H) ~ 8.5. 
The most recent spectroscopic observations of a sample 
of Hn regions i n M31 date from the late 1 990's and 
were presented by Galarza, Waltcr bos. fc Braunl (1 19991 ) and 
iBresoIin. Kennicutt. fc Garnettl (|l999h . The former use the 
R23 parameter and find a somewhat steeper oxygen gradi- 
ent, with a slope of — 0.06±0.03 dex kpc -1 . All these authors 
derive the chemical abundance gradient in M31 by applying 
empirical calibrations of bright-line ratios, as the tempera- 
ture sensitive auroral lines used to constrain electron tem- 
perature and emissivities of abundance-sensitive forbidden 
lines (i.e. through the direct method) remained undetected. 
In fact, it is somewhat surprising that, to-date, there are no 
auroral line detections for the determination of Hn region 
abundances through the direct metho d in M31, apart from a 
single Hn region (K932) observed bv lEsteban et all (|2009h . 

Other attempts to measure the oxygen abun- 
dance gradient in M31 come from the reanalysis 
of the datasets mentioned above (mainly those of 
the 1980's) and the application of different s tronq- 



Table 1. Journal of observations. 



line methods llZaritskv. Kennicutt. fc Huchra 
Vila-Costas fc Edmunds! Il992i ; Smartt et al 



1994 



2001 



Trundle et al.l [2002). The M31 oxygen abundance gra- 



dient slope derived by the different authors ranges from 
—0.013 dex kpc -1 to —0.06 dex kpc -1 , i.e. they differ by as 
much as a factor of four, mainly due to uncertainties and 
differences in the strong-line methods used. The central 



Mask 


AA 


Spectral range a 


Exposure Time 




(A, FWHM) 


(A) 


00 


1 


5.6 


3450-6000 


3x1800 + 2x1200 


1 


4.5 


5120-6840 


3x1800 


1 


9.3 


6450-10000 


2x1200 


2 


5.6 


3450-6000 


3x1800 + 2x1000 


2 


4.5 


5120-6840 


3x1800 


2 


9.3 


6450-10000 


2x1000 


3 


5.6 


3450-6000 


2x1800 + 2x1000 


3 


4.5 


5120-6840 


2x1800 


3 


9.3 


6450-10000 


2x1000 



a Approximate spectral range for a slit in the centre of the detector. 

oxygen abundance that results from these studies ranges 
between 1 and 3.5 times the solar value. 

Motivated by the desirability of a modern re-evaluation 
of the present-day chemical abundance in M31, we secured 
new spectra of H 11 regions in the Andromeda galaxy, in or- 
der to obtain nebular chemical abundances based on the 
detection of the faint auroral lines used for electron tem- 
perature determinations. The main aim of the current pa- 
per is to measure the oxygen abundance gradient from Hn 
regions in M31 and to compare with othe r metallicity indi- 
cators, in particular young, massive sta rs (|Venn et al.ll2000l ; 
iTrundle et al.ll2002l; ISmartt et al1l200ll ) and planetary neb- 
ulae l|Kwitter et al.ll2012i '). 

This paper is organized as follows. The spectroscopic 
observations, the data reduction procedure and the method- 
ology employed for meassuring emission line fluxes are de- 
scribed in Sec. [2] In Sec. [3] we derive the physical properties 
(electron densities and temperatures) of the H n regions and 
from these, the ionic and total chemical abundances are ob- 
tained (Sec. |4|. The radial oxygen abundance gradient is 
then calculated (Sec. [5]l from our direct measurements and 
from empirical metallicity calibrations applied to the Hn 
regions of our sample and to a compilation of data from 
other authors. A comparison with other metallicity indica- 
tors is presented in Sec.[|J] Finally, in Sec.[7]with discuss and 
summarize our results. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The spectroscopic observations of H 11 regions in the disk of 
M31 were carried out on 2005 September 29 with the Low 
Resolution Imaging Spectrometer (LRIS, lOke et all Il995l ) 
at the 10m Keck I telescope operated by the W.M. Keck 
Observatory on the Mauna Kea summit. The sky was pho- 
tometric with seeing ^0.8" at the beginning of the night. 
Conditions slowly got worse, having seeing ~1.5" at the end 
of the night. 

LRIS was used in multi-slit mode. Both the blue and red 
channels of the spectrograph were used simultaneously with 
a dichroic beam splitter (with 50% transmittance at 5091 A) 
for two different setups. The first setup was used to cover 
the spectral range 3500 A to 6840 A for a slit at the cen- 
ter of the detector with spectral resolution and dispersion 
~5.6A (FWHM) and 0.62A/pix, respectively, in the blue 
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Figure 1. Identification of the targets included in our LRIS Keck multi-slit masks plotted over an Ho continuum-subtracted image 
obtained at the 2.5m Isaac Newton Telescope at la Palma. The circles mark the brightest Ha knots along each slit, where the apertures 
for the one dimensional spectrum extractions were centred. 



Table 2. M31 adopted parameters. 



Parameter 


Value 


References 


RA (J2000) 


OOh 42m 44.33s 


Skrutskic ct al. (20061 


Dec (J2000) 


+41° 16' 07.50" 


Skrutskie et al. (20061 


Inclination 


77° 


Corbelli et al. (2010) 


Position angle 


38° 


Corbelli et al. (2010) 


Distance 


(744 ± 33) kpc 


Vilardell et al. (20101 


R25 (arcmin) 


(95 + 2) 


dc Vaucoulcurs ct al. (1991) 


R25 (kpc) 


(21 + 1) 





channel (with a 600 lines mm 1 grism blazed at 4000 A) 
and ~4.5A (FWHM) and 0.62A/pix in the red channel 
(with a 900 lines mm -1 grating blazed at 5500 A). The sec- 
ond setup was identical in the blue channel, but included 
a 400 lines mm -1 grating blazed at 8500 A in the red arm, 
covering roughly from 6450 A to 10000 A with 9.3 A FWHM 
resolution. The blue channel spectra were imaged on a mo- 
saic of two 2Kx4K Marconi (E2V) CCDs giving a plate 
scale of 0.135" pix -1 , the red channel spectra on a back- 
illuminated Tektronix 2Kx2K CCD, with 0.215" pix" 1 . 

Three slit masks were defined and used for the multi- 
object spectroscopy, centered approximately 17', 42' and 66' 
from the nucleus of M31. The masks cover 4.8' x 7.3' on the 
sky, and contained 8, 15 and 13 1.2"-wide slits, respectively. 



The slit lengths varied between 7" and 1.5'. The observa- 
tions were carried out ensuring that the slit position angle 
in each mask was within 20° of the parallactic angle for air- 
masses greater than 1.3, in order to minimize light losses 
due to differential atmospheric refraction. Table [T] contains 
a journal of the observations. 

The target Hn regions were selected from Ha 
continuum-subtracted images (see Fig. [l}, previously ob- 
tained with the Wide Field Camera of the 2.5m Isaac New- 
ton Telescope (INT) at La Palma (Spain) withi n the Local 
Group Census project (e.g. iMaerini et alJliooj ). An astro- 
metric calibration of these images with an accuracy of 0.3 
arcsec was se cured using stellar positions from the USN02 
l|Monet|[l998l ) catalogue. The positions of the selected Hn 
regions across the M31 disc are shown in Fig.[2]and their co- 
ordinate s are given i n Tab le EH Four targets are in common 
with the lBlair et"aH (|l982h sample (BA423, BA379, BA416, 
BA370; the last two were classified as SNRs by these au- 
thors). 



1 The Local Group Census is a narrow-band imaging 
survey of all the galaxies of the Local Group above 
Dec= — 30°(http:/ / w ww.ing.iac.es / ~rcorradi /LGC / \ . 
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Figure 2. Image of M31 showin g the location of the H II studied in this w ork (white circles ). Hn regions studied by other authors 
are indicated with p ink triangles l|Dennefeld &, K untb||l9 8l | ). violet diamonds l|Blair et el1.II 1982^) . yellow squa res (showing slit centers by 
iGalarza et al.|[l999] ). inverted orange triangles lBresolin et al.lll999h and cyan pentagon l lEsteban et alj|2009h . The image is a composite 
of GALEX (FUV: blue, NUV: green) and Spitzer (24 fj,m; red). 



2.2 Data reduction 

The data reduction was carried out using standard IRAF0 
techniques. After removing the bias level, a correction was 
applied to account for the different gains (e~/ADU) of the 
amplifiers of the red and blue CCDs. Afterwards the spec- 
trum corresponding to each slit was identified and sepa- 
rated out from the 2-dimensional slit mask spectra of sci- 
ence and calibration frames, and in subsequent tasks the 

2 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 



spectrum from each slit was treated as a separate spec- 
trum. The data reduction was then carried out using stan- 
dard procedures for lo ng-slit spectra: natfi elding, cosmic-ray 
removal (lacos_spec, Ivan Dokkum|[200lh . removal of geo- 
metrical distortions in the spectral and spatial directions, 
wavelength calibration, extraction of one dimensional spec- 
tra within defined apertures and combination of the individ- 
ual spectra. 

The final spectra were flux calibrated using observa- 
tions of the spectrophotometric standards BD+28 4211, 
G191-B2B, Hiltner 600 and GD71 that were obtained at 
the beginning, middle and end of the night with a 8 .7" - wide 
slit. T h e standard star fluxes have been obtai ned from | Oke 
l|l990l) . iBohlin. Colina. fc Finlevl (|l995l ). and lHamuv et a~ 
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Table 3. Hll region equatorial coordinates and galactocentric distances. 



ID^ 


RA 


Dec 


Galactocentric distance'' 


Other ID£ 




(J2000) 


(kpc) 




1 


00:41:43.67 


+41:08:40.27 


5.05 




2* 


00:41:19.94 


+41:08:21.38 


8.16 


BA 419, P 318 


3t 


00:41:13.96 


+41:07:35.09 


8.64 


BA 422, P 310 


4 


00:41:14.54 


+41:07:17.25 


8.42 


P 308 




00:41:35.65 


+41:06:53.86 


5.42 


BA 416 (SNR candidate), P 325* 


6 


00:41:47.42 


+41:06:39.13 


3.93 


BA 415 (SNR candidate), P 327 


7 


00-41 -34 30 


+41:05:52.14 


5.23 


P 317 


at.* 


00:41:30.17 


+41:05:01.51 


5.47 


BA 423, P 309* 


9 


00:40:26.10 


+40:45:08.45 


8.85 


P 182 


10 


00:40:40.52 


+40:44:26.75 


8.58 




11 


00-40-27 26 


-I-40-44-00 71 

I T:U . t:t: * UU . i _1_ 


8.93 


BA 464, P 176 


12 


00:40:26.08 


+40:43:43.06 


9.01 


P 173 


13 


00:40:33.81 


+40:43:33.39 


8.86 


P 179 


14* 


00:40:28.77 


+40:43:11.74 


9.02 


P 171 


15 


00:40:33.35 


+40:43:03.86 


8.97 


P 174? 


16 


00:40:32.07 


+40:42:50.72 


9.04 


P 174? 


17* 


00-40-30 81 


4-40-42-30 19 


9.13 


BA 310, P 167 


18 


00:40:29.81 


+40:42:16.84 


9.19 


— 


19* 


00-40-31 55 


+40:41:21.40 


9.41 


BA 312, P 161 


20 


00:40:31.11 


+40:39:52.24 


9.85 




21 


00:40:32.94 


+40:39:18.35 


10.07 


BA 318, P 150 


22* 


00:40:30.48 


+40:39:09.33 


10.07 


P 145 


23 


00:40:32.57 


+40:38:57.02 


10.18 


BA 319, P 149 


24 


00:39:18.84 


+40:22:26.04 


14.53 


BA 375, P 15 


25* 


00:39:18.77 


+40:21:58.66 


14.66 


BA 379, P 7* 


26* 


00:39:18.77 


+40:21:58.66 


14.66 


BA 379, P 7* 


27 


00:39:21.79 


+40:21:44.31 


14.74 


BA 378, P 9 


28 


00:39:41.55 


+40:21:33.54 


15.25 


BA 358, P 32 


29*'* 


00:39:31.53 


+40:21:15.81 


15.06 


BA 368, P 21 


30 


00:39:47.97 


+40:21:11.07 


15.67 


BA 357, P 34 


31 


00:39:30.12 


+40:21:06.05 


15.09 


BA 373 


32*. * 


00:39:44.93 


+40:20:43.43 


15.73 


BA 360 


33* 


00:39:42.92 


+40:20:39.90 


15.67 


BA 362 


34* 


00:39:33.59 


+40:20:26.06 


15.42 


BA 374, P 17 


35* 


00:39:37.44 


+40:20:11.49 


15.65 


BA 371, P 20 


36 


00:39:35.28 


+40:19:43.06 


15.76 


BA 377, P 11 


37 


00:39:39.13 


+40:19:21.71 


16.06 


BA 376, P 14 



a Nebulae with auroral line detections are indicated with *. A double dagger (*) indicates supernova remnant candidates, 
k Galactocentric distances c omputed from the M 31 orientation p arameters shown in Table [2] 

c IDs starti ng with BA fromlBaade fc Ard lll964l) and with P from lPellet et al.l ( Il978ll . Nebulae marked with * have been spectroscopically 
observed bv lBIair etHl l|l982 V 

{ 19941 ). A standard mean optical/near-IR extinction curve 
for Mauna Kea was adopted for the atmospheric extinction 
correction (from the UKIRT web page^f]). 

Given the limited near-IR range of most of 
the spectrophotometric standard star spectra, only 
data below 9200 A are useful. The [Sm]A9532 fluxes 
have been calculated from the theoretical ratio 
[S m]A9532/[S m]A9069=2.44. 



3 http://www.jach.hawaii.edu/UKIRT/astronomy/utils/exts.html 
This extinction curve comes from CFHT Bulletin 
llBeland. Boulade, k. Davidgel ll98Sl) and the CFHT observer's 
manual. 



2.3 Emission line flux measurements and 
interstellar extinction correction 

The emission line fluxes of the Hn regions were measured 
with the SPL0T 1RAF task, by integrating all the flux in 
the emission line between two given limits over the local 
continuum level defined interactively. 

We analysed the spectra of 37 nebulae. Two of them 
correspond to different extractions of the same slit, in which 
two bright knots of ionized gas emission (objects 25 and 26 
in Table|3]) were identified. For 13 targets we were able to de- 
tect one or several of the following auroral lines, [O iii]A4363, 
[Nll]A5755, [Sm]A6312, and [Oll]A7325, which allow the de- 
termination of the electron temperature of the different ions, 
although four of the targets turned out to be SNR, as de- 
scribed below. The other 24 targets have lower signal-to- 
noise spectra and no auroral lines were detected. We iden- 
tified two SNR candidates within this subsample. For the 
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remaining 22 Hn regions we concentrate here only on the 
brightest emission lines, which will be used for the determi- 
nation of the oxygen abundances through strong-line meth- 
ods (Sect. EHj) . 

The three different spectral ranges (blue, red and near- 
infrared) were scaled in order to bring them to the same 
relative flux calibration. Several emission lines in common 
between the red and the near-infrared spectral ranges (typ- 
ically Ha, [Nn]AA6548,6583, Hel A6678, [S n]AA6717, 6731) 
provided the scaling factor. There are no lines in common 
between the blue and the red spectral ranges. In this case 
a preliminary scaling factor was obtained from the con- 
tinuum level of both spectra. The final scaling factor was 
late r determined by imposin g the expected Ha/H/3 flux ra- 
tio (jStorev fc Hummed Il995h for the measured or assumed 
electron temperature, after the interstellar extinction and 
stellar absorption corrections had been applied (see below). 

The interstellar reddening was determined from the 
Balmer series hydrogen recombination lines flux decrement. 
As the fluxes of nebular Balmer emission lines are affected 
by absorption by the underlying stellar population, a correc- 
tion needs to be applied in order to obtain the interstellar 
extinction to properly correct our measurements. We have 
assumed that the equivalent width of the stellar absorption 
(EWjij) is the same for all Balmer lines. Then, an iterative 
procedure was performed: the extinction coefficient c(H/3) 
was determined from the comparison of the measured fluxes 
of the H7 and H<5 lines relative to H/3 and the expected ratios 
from t he case B recombination given bv lStorev fc Hummer! 
1 19951 ) for an elect ron temp e ratur e of 7500 K and density 
of 100 cm -3 . The iHowa rth ( 1983) param etrization of the 
interstellar reddening law bv lSeatonl l| 19791 ) was adopted for 
R^=3.1. The EWais was then adjusted in successive itera- 
tions until the c(H/3) value obtained from the H7/H/3 and 
H5/H/3 line ratios were in agreement. After the determina- 
tion of the electron temperature (see Sect. [3} the iterative 
process was repeated. We checked that higher order recom- 
bination lines of the Balmer series measured in our spectra 
were in agreement (within errors) with the theoretical pre- 
dictions. The extinction coefficient c(H/3) obtained varies 
between 0.24 and 1.0 in our Hn region sample, while the 
EWats was found to be in the range 1.5-3.2 A for the high- 
est S/N spectra, where we detected the auroral lines. For 
the lowest S/N spectra, the extinction coefficient c(H/3) was 
typically determined from the H7/H/3 ratio alone (and in 
this case the EW„h was assumed to be zero) or from both 
Hy/H/3 and H5/H/3 when possible. c(H/3) for these regions 
ranges from 0.0 to 1.03. 

The final values of c(H/3), EW a j s , and the blue-red scal- 
ing factor were further checked in the highest S/N Hn re- 
gion spectra by compar ing the observed He I line ratios wit h 
theoretical predictions (iBeniamin. Skillman. fc Smitslll999h . 
The observed and theoretical Hel line ratios A5876/A4471, 
A6678/A4471 and A5876/A4922 differ by less than 15% in all 
targets, except for #5, for which differences are larger (of 
order —40%) . We note that the observational errors for these 
ratios are also higher (—30%) for this target, which turned 
out to be a SNR (see Sect. [2^4)1 . The same check was done 
with the Balmer ratio He/H/3 and Paschen series line ratios 
to H/3 when available (using mainly Pa 9, Pa 10, Pa 11, Pa 
12) . In this case the observed and theoretical line ratios differ 
by less than 20% (typically -5-10%), except for the He/H/3 



ratio in targets #5, #25, #26 and #35 for which the dif- 
ference observed-theoretical is larger, —25-30%, except for 
target #5 (SNR) in which He is twice the value expected 
for an Hn region possibly due to blendi ng with enhanced 
[Nem] jStupar. Parker, fc Filipovicj|2008r ). 

Table 2] shows the line fluxes corrected for reddening 
and normalized to H/3 for the 13 nebulae where auroral lines 
fluxes have been measured. For the remaining targets (with 
no detections of auroral lines) we give in Table [S]the bright- 
est line fluxes corrected for reddening and normalized to H/3. 
The coordinates and galactocentric distances for both sets 
of regions are given in Table [3] 

The uncertainties quoted in Tables [4] and [5] include 
(added in quadrature): statistical errors measured with 
SPL0T, flatfielding errors, the uncertainty in the flux cal- 
ibration, and the uncertainty in the determination of the 
flux scaling factors between the different spectral ranges and 
in the determination of the extinction coefficient. Typically, 
errors in the line flux ratios are of order —4-8% for lines 
brighter than Fhp/2, and —6-16% for fainter lines, for the 
Hn regions with detections of auroral lines. For the regions 
with lower S/N the errors in the line flux ratios relevant for 
strong-line methods for abundance determinations are of or- 
der -5-20% for [Oll]/H/3, -5-40% for [O m]A4959/H/3, -4- 
30% for [Om]A5007/H/3, and -5-20% for [N ll]A6583/H/3. 

Four of our t argets are located within —3" of the posi- 
tions reported by iBlair et alj (| 19821 . see Table |3J . We have 
compared the fluxes published by these authors with the 
fluxes of our brightest lines. The measurements agree within 
15% for lines brighter than Fj^/3 across the whole spectral 
range m common (from 3700 A to 6730 A). The differences 
are higher for [S ll]AA6717, 6731 in #5 (BA416), -35%, and 
for [O n]A3727 in #25 (BA379), -50%. For line ratios fainter 
than Fhp/3 the differences increase up to —60% (in both di- 
rections), with no dependence on wavelength. We note that 
IBlair et al.l (| 19821 ) state that their error in faint line fluxes 
(<0.5F_h^) may be as large as 50%. 



2.4 Supernova remnants 

Targets # 5 (BA416) and # 29 (BA370) are classified as 
SNRs by IBlair et al.l (| 19821 ). For both nebulae we find 
[Sn]/Ha>0.5, indicating the presence of shocks. Our mea- 
sured line ratios for BA416 position this region in the area 
of the Hq/[Sii] v s. Ha/ [Nil] diagnostic diagram that cor- 
responds to SNR (|Sabbadin. Minello. fc Bianchinill 19771) . in 



agreement with the conclusion by Bl air et al. (|l982h . Ac- 



cording to our measurements, the location of BA370 in the 
same diagnostic diagram does not allow us to clearly estab- 
lish the nature of this nebula. However, we measure a sub- 
stantial [Oi]AA6300,6360 flux ( [ O i] /R/3 > 0.4) . Following 
the lin e ratio diagnostic criteria of lFesen. Blair, fc Kirshner] 
jl985l ) ([On] /H/3 vs. [Oi]/H/3) both BA416 and BA370 
are identifi e d as SNRs. We note that, as mentioned above, 
IBlair et al.l i|l982i ) reported a -35% higher [S n]AA6717, 6731 
flux for BA416 than determined from our observations. 
Given that the center position of our slit does not perfectly 
coincide with theirs, we cannot exclude some dilution of the 
SNR emission by surrounding photoionized gas emission in 
our spectrum. 



Table 4: Reddening-corrected line fluxes in nebulae with auroral line 
detections. 
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Table 4: Continued. 
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11±2 


19±1 


23±3 


19.3±1.3 


21±2 


22±2 


8.00±2 


21±2 


16±2 


20±3 


13.4±1.2 




9229 


2.5±0.4 






2.43±0.13 


2.9±0.4 




2.4±0.2 


2.0±0.2 




2.3±0.4 


2.0±0.8 


2±2 




m 




0.21±0.04 


0.31±0.02 


0.76±0.1 


0.67±0.02 


0.55±0.10 


0.44±0.02 


0.43±0.05 


0.42±0.04 


0.34±0.04 


1.00±0.05 


0.39±0.05 


0.54±0.04 


0.44±0.02 






189.6±0.4 


198.7±0.4 


166.6±1.1 


156.90±0.14 


226.1±0.4 


74.9±0.3 


123.4±0.1 


377.6±0.6 


174.1±0.8 


908±2 


186.0±0.2 


82.3±0.2 


184.0±0.5 


ipy* 




11.1±0.5 


13.0±0.5 


3.9±0.2 


142±5 


21±1 


9.8±0.3 


53±2 


25.5±1.4 


3.5±0.1 


26±1 


35.7±1.2 


17.3±0.8 


9.4±0.3 



dicates SNRs candidates (see Sect. [2~3t for details. 

'(H/3) is the total extinction— corrected H/3 flux in units of 10 — 15 erg cm -2 s _1 . 
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Figure 3. Comparison between the electron temper- 
ature obtained from the [O III] A4363/[0 III] AA4959, 5007 
line ratio and the electron temperatures ob- 
tained from [Siii]A6312/[Siii]AA9069,9532 (top), 
[O n] A7325/ [O il] AA3727.3729 (middle) , and 

[Nii]A5755/[Nii]AA6548,6584 (bottom). The cy an circle cor- 
respon ds to the Hn region K932, studied by lEsteban et al.l 
(2009). The green straight lines indicate the relat ionship between 
different temperatures predicted by the iGarnetd jl992h models. 



Figure 4. Comparison between the electron temper- 
ature obtained from the [N II] A5755/[N II] AA6548,6584 
line ratio and the electron temperatures ob- 
tained from [Siii]A6312/[Sni]AA9069,9532 (top), 
[On]A7325/[On]AA3727,3729 (middle) and 

[Sii]A4072/[Sii]AA6717,6731 (bottom). The cy an circle cor- 
respon ds to the H II region K932, studied by lEsteban et al.l 
[120091) . 



Four additional targets (#8, #14, #22 and #32) are 
SNR candidates. Targets #8 and #32 have a [Sll]/Ha ra- 
tio consistent with an Hn region, but their spectra present 
numerous metal lines in the blue range (mainly from Fe 
and Mg), as well as relatively strong [O i] AA6300, 6360. In 
addition, their [Om]A4363 emission appears to be contam- 
inated by [Feii]A4359. We suspect that these targets are 
young SNRs. 

For targets with lower signal-to-noise spectra, we have 
tentatively identified as SNRs those with [Sll]/Ha>0.4 
and substantial ]OilA63 00 emission, [O i]A6300 > 0.02, 
l|Dennefeld fe Kunthlll98"lT ). According to this criteria, tar- 
gets #14 and #22 are SNR candidates as well. 

The six targets discussed in this section are removed 
from our sample and will not be considered in the remainder 
of the paper, since the application of Hn region metallicity 
diagnostics would not be appropriate. For completeness we 
report their line fluxes in Table 0] (#5, #8, #29 and #32) 
and in Table [5] (#14 and #22). 



3 PHYSICAL CONDITIONS OF THE IONIZED 

GAS 

The physical conditions of the ionized gas (electron tem- 
perature and density) have been determined from colli- 
sionally excited line r atios. The temden IR AF task within 
the NEBULAR package l|Shaw fe Dufourl ll995) was employed. 



The temden task is an extension o f the FIVEL program 
|De Robertis. Dufour. fe Hund[l987l) for a larger set of ions, 
emission lines and number of levels (between 5 and 8 depend- 
ing on the ion). The atomic parameters (transition prob- 
abilities and collisional strengths) for the ions of interest 
were updated with mo re recent data as shown in Table 5 of 
iBresolin et all [|2009ah . 

The derivation of the electron temperature of an H II re- 
gion requires the detection of auroral lines which arise from 
the transition from the second lowest excited level to the 
lowest excited level. Auroral lines are generally faint and 
difficult to detect, with strongly decreasing intensities with 
increasing gas metallicity. This explains why, to our knowl- 
edge, there is surprisingly only one previous published mea- 
surement of electron temp erature for an H H reg ion in M31: 
region K932, analysed by lEsteban" et all lj2009l l from Keck 
HIRES data. For nine H II regions of our sample, one or more 
of the auroral lines [Om]A4363, [Sm]A6312, [Oll]A7325, 
[Nn]A5755 and [Sii]A4072 have been detected. The ratio 
of these lines to stronger nebular lines of the same ions are 
highly dependent on electron temperature and can thus be 
employed for measuring it. The auroral line [Nll]A5755 was 
detected in all nine Hn regions, while [Sm]A6312 was mea- 
sured in only five. This is in part due to the fact that this 
emission line, due to the radial velocity of M31, is within 
2 A of the telluric line [O i] A6300 A, and its firm detection is 
therefore subject to higher uncertainties. 

Density and temperature of the ionized gas have been 
simultaneously obtained from an iterative process. The tern- 
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Table 5. Reddening-corrected line fluxes relative to H/3 (with H/3=100) for the Hll regions with no detections of auroral lines. 



ID [On] [Oni] [Om] [Nn] Ha [Nn] [Sn] [Sn] F(H/3) c(H/3) 
3727 4959 5007 6548 6562 6583 6717 6731 (1(T 15 erg s^cm -2 ) (mag) 



1 








44 ±2 


290 ± 20 


149 ±8 






0.20 ±0.08 


0.00 


i 


140 ±30 






20 ±4 


290 ± 60 


60 ±13 


27 ± 6 


19 ±4 


0.5 ±0.1 


0.64 


6 


122 ±6 


4.7±0.5 


11.6±0.7 


33 ±2 


300 ± 20 


122 ±7 


40 ±2 


25 ±2 


1.43 ±0.05 


0.00 


7 


230 ± 20 


2.1±0.9 


3.5±1.0 


51 ±5 


290 ± 30 


161 ± 15 






1.09 ±0.04 


0.45 


9 


200 ± 20 


12 ±2 


35 ±3 


38 ±4 


290 ± 30 


123 ±12 


77 ±8 


53 ±5 


0.42 ±0.02 


0.00 


10 


150 ±20 




7±2 


45 ±6 


290 ± 40 


130 ± 16 


48 ±6 


33 ±4 


1.00 ±0.04 


0.58 


11 


150 ±10 


10.3 ±0.6 


29.5 ±1.4 


31 ±2 


290 ± 20 


93 ±7 


30 ±2 


21 ±2 


3.9 ±0.1 


0.19 


12 


149 ± 12 


6.2± 1.4 


18 ±2 


40 ±4 


290 ± 30 


120± 11 


72 ±7 


46 ±5 


0.232 ±0.009 


0.00 


13 


97 ±8 


6.5±0.8 


21 ±1 


26 ±2 


290 ± 30 


75 ±7 






1.15 ±0.04 


0.36 


14* 


390 ± 50 




16 ±3 


29 ±6 


280 ± 50 


100 ±20 


90 ±20 


69 ±13 


0.56 ±0.04 


0.59 


15 


180 ±20 


12 ±4 


27 ±4 


30 ±5 


280 ± 50 


80 ±13 


47±8 


29 ±5 


0.68 ±0.03 


0.64 


16 


94 ±8 




4.6 ±0.8 


18 ±2 


290 ± 30 


67 ±6 


36 ±4 


24 ±2 


0.58 ±0.02 


0.00 


18 


193 ±20 






39 ±4 


300 ± 30 


124 ±12 


72 ±7 


50 ±5 


0.68 ±0.03 


0.68 


20 


130 ±8 


6.6±0.6 


18 ±1 


25 ±2 


290 ± 20 


73 ±5 


31 ±2 


21 ±2 


8.9 ±0.3 


0.90 


21 


126 ±7 


3.2 ±0.4 


7.9±0.5 


27 ±2 


290 ± 20 


83 ±5 


29 ±2 


20.2 ±1.4 


8.5 ±0.3 


0.98 


22* 


220 ±20 


23 ±2 


62 ±3 


38 ±2 


290 ± 20 


123 ±8 


64 ±4 


46 ±3 


2.03 ±0.07 


0.61 


23 


130 ±8 


6.6±0.6 


18.0 ±1.1 


25 ±2 


290 ± 20 


73 ±5 


31 ±2 


21 ±2 


8.9 ±0.3 


0.90 


24 


300 ± 20 


22.1±1.2 


64 ±3 


31 ±2 


290 ± 20 


91 ±6 


48 ±3 


33 ±2 


2.80 ±0.09 


0.42 


27 


210 ±20 






19 ±2 


280 ± 30 


68 ±6 


42 ±4 


27 ±3 


3.8 ±0.1 


0.79 


28 




74.6 ±3.6 


220 ±10 


13.3 ±0.9 


290 ± 20 


43 ±3 


15±1 


10.2±0.7 


14.4 ±0.05 


1.03 


30 






7.3 ±0.3 


20.6 ±1.3 


290 ± 20 


62 ±4 


25 ±2 


16.9 ±1.1 


2.52 ±0.08 


0.19 


31 


290 ± 20 


11.2±0.8 


35 ±2 


28 ±2 


290 ± 20 


83 ±5 


38 ±3 


28 ±2 


5.9 ±0.2 


0.57 


36 


253 ±19 


6.3±0.3 


15.6 ±0.7 


21 ±2 


280 ± 30 


56.5 ±5 






7.7 ±0.3 


1.03 


37 


257 ±13 


16.2 ±0.8 


48 ±2 


25 ±2 


280 ± 20 


75 ±4 






13.3 ±0.5 


0.74 



* SNRs candidates (see Sect. 12.41 1 for details. 



Table 6. Measured Hll region electron temperatures (in K) and densities (in cm -3 ). The 
auroral line intensities [Oll]A7325 and [Nn]A5755 have been corrected for recombination 
contamination following iLiu et al ] l|2000l l. 



ID 


T[0 in] 


T[S in]* 


T[On] 


T[Nn] 


T[S n] 


Ne 




4363 


6312 


7325 


5755 


4072 


[S n]AA6717/6731 




4959+5007 


9069+9532 


3727 


6548+6583 


6717+6731 


2 




7100 ± 300 


7200 ± 500 


7500 ± 400 


7200 ± 400 


66 


3 








8300 ± 300 




19 


17 


8400 ± 200 


8300 ± 500 


9300 ± 700 


7800 ±400 


7500 ± 500 


38 


19 








7200 ±400 


6900 ± 200 


24** 


25 


9300 ± 200 


9300 ±400 


10600 ± 1700 


9600 ± 400 


7800 ± 300 


28 


26 


9500 ± 300 


9800 ±400 


11200 ±800 


9900 ±400 


7300 ±400 




33 






13300 ± 1000 


8800 ± 300 




24 


34 






14400 ± 900 


10800 ± 500 




14** 


35 


8800 ± 200 


10100 ±600 




9800 ± 800 


7300 ± 300 


80 



* Assuming [S m]A9532/[S in] A9069=2.44. 

** [Sn]6717/[Sn]6731 line ratios slightly above 1.41 (1.44, 1.43 and 1.45 respectively); it 
has been assumed that [S n] A6717/[S n] A6731=1.41. 



perature was initially set to T e =7500 K and a prelimi- 
nary value for the electron density N e was derived from 
the [Sii]A6717/[Sii]A6731 line ratio. This initial value of N e 
was then used for computing T e for the different auroral-to- 
nebular line ratios. The process was repeated until conver- 
gence. An initial estimate of the ionic abundances was then 
obtained (as described in Sect. I4.l( l and from these we esti- 
mated the contamination of the auroral lin es [Nn]A5755 and 
[On]A7325 due to recombination following ILiu et ail (|200fJ ) 
(their Eqs. 1 and 2) and assuming N++/N+ ~ 0++/0+ to 
obtain N ++ /H + . The electron temperatures and densities 
were recomputed until convergence after the correction for 
recombination contamination was applied. This correction 
is small and varies between 0.1 and 0.9% for [N n]A5755/H/3 
and between 0.3 and 3.7% for [O n] A7325/H/3. We note that 



in a number o f regio ns T[0 ll] is above the range of validity of 
the lLiu et all (|2000l l formula for [On]A7325/H^ (5000 K < 
T e < 10000 K). In any case, the errors in the derived T[0 n] 
value for these Hn regions are significantly larger than the 
corresponding correction. 

The electron temperatures are shown in columns 2 to 6 
of Table [6] The errors reported for T e arise from error prop- 
agation of the relevant line flux errors and are in the range 
200-600 K for all ions, except in the case of T[0 ll] , for which 
the errors are larger (400-1700 K). This is mainly due to the 
observational uncertainty that involves the measurement of 
the [O n] multiplet at 7325 A, located in a spectral region of 
strong OH airglow emission. 

Figs. [3] and [4] show a comparison between the different 
ionic temperatures we measured and the results one would 
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obtain from model predictions bv lGarnettl(|l992h , frequently 
used for predicting ionic temperatures when only one or two 
auroral lines are detected. Fig.[3]shows the relations between 
T[Sm], T[Oll] and T[Nn] with respect to T[Om]. For most 
of the regions with a T[0 ill] determination, the iGarnettl 
(1992) temperature relations 



Table 7. Adopted electron temperatures for the two— zone scheme 
of the H II region ionization structure. 



T[Sm] = 0.83T[Oin] + 1700K 

T[0n] = T[Nn] = 0.70T[Oin] + 3000K 



(1) 
(2) 



appear to be satisfactorily followed by the empirical 
data, especially for the case of T[Sm] vs. T[Oiii]. 
In the other two cases the data points deviate more 
significantly from the theoretical relations. This was 
also found by oth er authors with a larger sample of 
Hii regions (e.g. Kennicutt. Bresolin. fc Garnett! 120031 : 



iPerez-Montero fc Dfadl2003| 1 

Fig. U shows the relations between the different ionic 
temperatures with respect to T[Nn], and highlights a resul t 
already found by other authors (e.g. iBresolin et al.l [2005h : 
T[Oll] seems to overestimate the temperature with respect 
to T[Nll], while the opposite happens to T[Sn], which gives 
always lower temperatures than T[Nn]. However, T[Sm] 
seems to be well correla ted with TJNii], following closely 
the model predictions bv lGarnettl (|l992h . 

The electron densities N e are shown in column 7 of Ta- 
ble [6] All Hii regions are in the low-density regime, with 
densities below 100 cm -3 . For three Hii regions we mea- 
sured a line ratio [S n] A6717/[S Il]A6731 slightly above the 
theoretical limit. For those regions we assumed a line ratio 
equal to 1.41 (equivalent to the theoretical limit). 



4 CHEMICAL ABUNDANCES 
4.1 Ionic abundances 

The calculation of ionic abundances from collisionally- 
excited lines requires a good knowledge of the electron tem- 
perature. An onion-model is commonly adopted to describe 
the ionization structure of an H II region, with a number of 
zones of different T e where atomic species of similar excita- 
tion coexist. In light of the correlations found in Sec. [3] be- 
tween the T e obtained from different atomic species, we have 
adopted a two-zone representation for describing the tem- 
perature stratification of the Hii regions: a low-excitation 
zone, characterized by T[Nll], where ions of low ionization 
potential coexist (0 + , S + , N + ) and a high-excitation zone 
for the ions with higher ionization potential (0 ++ , S ++ , 
Ar ++ , Ne ++ ), characterized by T[Om]. For the Hii regions 
for which we could not measure T[Om], we have used the 
average of our measured value for T[Siii], and the value 
for TjOiiil obtained from the measured T[Nll] by invert- 
ing Eq.[2l (|Garnettlll992l ). When T[Sm] was not available, 
T[0 ill] has been obtained from T[N n] and the mentioned re- 
lation (Eq.[2]). TableQcontains the temperatures adopted for 
the ionic abundance calculation for each H II region. The er- 
rors quoted in the table correspond to the errors in the mea- 
sured values for T[Nn] and T[Om], as described in Sect.|3j 
when these temperatures were available. When the Garnett 
relations were used, the uncertainty in T[Om] comes from 
error propagation of the uncertainty in T[Nll], T[S III] and 



ID 


T[0+, N+, S+] 


T[S++, Ar++, 0++, Ne++] 




(K) 


(K) 


2 


7500 ± 400 


6800 ± 400 


3 


8300 ± 300 


7600 ± 600 


17 


7800 ± 400 


8400 ± 200 


19 


7200 ± 400 


6000 ± 700 


25 


9600 ± 400 


9300 ± 200 


26 


9900 ± 400 


9500 ± 300 


33 


8800 ± 300 


8300 ± 600 


34 


10800 ± 500 


11100 ±800 


35 


9800 ± 800 


8800 ± 200 



an additional term added in quadrature, which is an esti- 
mate of the uncertainty in the scaling relation (±400K for 
T[N n] vs. T[0 m]) obtained from wo rks with a large number 
of Hii regions (Bres olin et alj|2009al ). 

The ionic IRAF task within the NEBULAR package 
dShaw fc Dufourlll995l ) was employed for obtaining the ionic 
abundances of 0+, ++ , N+, S + , S ++ , Ar ++ and Ne ++ . 
They were estimated from the relevant reddening-corrected 
emission line ratios from those species (Tabled and the cor- 
responding electron temperatures of the adopted two-zone 
scheme (Table • Table [5] shows the resulting ionic abun- 
dances with respect to H + . 

4.2 Total abundances 

The total abundance relative to hydrogen of a given ele- 
ment is given by the sum of the abundances of all its ions. 
However, not all ions have emission lines in our spectra 
and we have to correct for these unseen ionization stages 
using ionization correction factors (ICFs). We have used 
the following common assumptions: 

Oxygen. Hell A4686 is only detected in one target 
(#8) which is a SNR. This implies that all H II regions are 
of sufficiently low excitation that the amount of oxygen in 
ionization stages higher than ++ is negligible. We have 
therefore set O/H = 0+/H+ + 0++/H+. 

Nitrogen and Neon. For these elements we have as- 
sumed that N/0 = N + /0 + and Ne/O = Ne ++ /0 ++ , which 
result from i onization potential c onsiderations of the ions 
involved (see IBresolin et al.ll2009al ). 

Argon. Only Ar ++ is observed in our optical 
spectra. However, the amount of Ar + in low ex- 
citation Hu regions (our case) can be important 
l|Bresolin. Garnett. fc Kennicuttll2004l ). We have employed 
the ICFs obtained from photoionization models by Izo- 
tov et al. (2006). These authors proposed different pa- 
rameterizations depending on metallicity. Those for high 
(12+log(0/H) > 8.2) and intermediate (7.6 < 12+log(0/H) 
^ 8.2) metallicity cover our range of oxygen abundances: 



ICF(Ar++) = 



0.285?; + 0.833 + 0.051/?; 
0.517« + 0.763 + 0.042/?; 



intermediate Z 
high Z 



where v = 0+/(0++0 ++ ). 



Sulfur. In the case of low excitation H II regions the 
ICF to correct (S + + S ++ )/H + for the presence of higher 
ionization sulfur stages is expected to be small. We have used 
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Table 8. + /0 and ionic abundances (columns 3 to 9) given as 12 + log(X/H), where X is the ion in the 
table header. 



ID 


o+/o 


o+ 


o++ 


N+ 


S+ 


S++ 


Nc++ 


Ar++ 


2 


0.8 ±0.3 


8.50±0.14 


8.02±0.13 


7.56 ±0.09 


6.30 ±0.08 


6.95 ±0.07 


7.21 ±0.16 


6.22 ±0.09 


3 


0.8 ±0.2 


8.16 ±0.07 


7.66 ±0.14 


7.41 ±0.04 


6.16 ±0.04 


6.78 ±0.08 




6.06 ±0.10 


17 


0.7 ±0.2 


8.25±0.11 


7.92 ±0.04 


7.40 ±0.07 


6.07 ±0.06 


6.73 ±0.02 


7.03 ±0.05 


6.11 ±0.03 


19 


0.7 ±0.3 


8.49 ±0.15 


8.06 ±0.31 


7.70 ±0.09 


6.46 ±0.08 


7.03 ±0.16 




6.36 ±0.21 


25 


0.4±0.1 


7.86 ±0.07 


8.08 ±0.03 


6.83 ±0.05 


5.68 ±0.04 


6.59 ±0.02 


7.34 ±0.04 


6.03 ±0.02 


26 


0.5 ±0.1 


7.95 ±0.08 


8.00 ±0.05 


6.98 ±0.05 


5.85 ±0.04 


6.58 ±0.03 


7.26 ±0.06 


6.01 ±0.03 


33 


0.8 ±0.2 


8.17±0.08 


7.65 ±0.14 


7.29 ±0.05 


6.18 ±0.05 


6.67 ±0.08 




5.91 ±0.09 


34 


0.6 ±0.1 


7.90 ±0.08 


7.68 ±0.10 


7.19 ±0.05 


6.13 ±0.05 


6.24 ±0.06 


6.63±0.12 


5.93 ±0.07 


35 


0.3 ±0.1 


7.80 ±0.16 


8.27 ±0.03 


6.84 ±0.09 


5.72 ±0.09 


6.64 ±0.02 


7.55 ±0.04 


6.13 ±0.02 



Table 9. Total abundances. 



ID 


12+log(0/H) 


log(N/0) 


log (S/O) 


log (Ar/O) 


log (Nc/O) 


Hc/H 


R/R.25* 


2 


8.62±0.11 


-0.95 ±0.17 


-1.58 ±0.13 


-2.32±0.14 


-0.81 ±0.20 


0.098 ±0.027 


0.40 


3 


8.28 ±0.06 


-0.76 ±0.08 


-1.41 ±0.09 


-2.14 ±0.12 




0.111 ±0.024 


0.42 


17 


8.42 ±0.08 


-0.86 ±0.13 


-1.59±0.08 


-2.23 ±0.08 


-0.89 ±0.06 


0.111 ±0.017 


0.44 


19 


8.63 ±0.14 


-0.79 ±0.17 


-1.49 ±0.19 


-2.19±0.25 




0.083 ±0.034 


0.46 


25 


8.28 ±0.04 


-1.02 ±0.09 


-1.58±0.04 


-2.23 ±0.04 


-0.74 ±0.05 


0.088 ±0.005 


0.71 


26 


8.28 ±0.04 


-0.98 ±0.09 


-1.59±0.05 


-2.23 ±0.06 


-0.73 ±0.07 


0.090 ±0.006 


0.71 


33 


8.29 ±0.07 


-0.89 ±0.10 


-1.49 ±0.09 


-2.30 ±0.12 




0.074 ±0.015 


0.76 


34 


8.10 ±0.06 


-0.71 ±0.10 


-1.60±0.07 


-2.14 ±0.09 


-1.05±0.16 


0.089 ±0.004 


0.75 


35 


8.40 ±0.05 


-0.96 ±0.18 


-1.60±0.05 


-2.24 ±0.05 


-0.72 ±0.05 


0.087 ±0.005 


0.76 



* Deprojected galactoccntric distance normalized to R.25- See Table [2] for the M31 adopted parameters. 



regions with considerable error bars, the curve with a = 2.5 
and log(S/0) = — 1.6 follows the data rather well. We have 
then kept this parameter choice, in a greement with previous 
work on extragalac tic Hn regions |Kennicutt et all 1 120031 ; 
iBresolin et al.ll2004h . This equation yields corrections of 16- 
30% for the highest excitation Hn regions of the sample 
(0+/0~ 0.25-0.4), 2-9% for Hn regions with O+/O-0.5- 
0.7, and corrections around 1% for the lowest excitation re- 
gions with O+/O~0.8. 

The total abundances of O, N, S, Ar and Ne, obtained 
as described above, are summarized in Table [5] 

4.3 Abundance trends with log(0/H) and 
galactocentric radius 

Fig. [6] shows the total abundances of nitrogen, sulfur, neon 
and argon with respect to oxygen as a function of the oxy- 
gen abundances. The N/O ratio does not show the typ- 
ical increase with 12+log( Q/H) expected for a pure sec - 
ondary origin of nitrogen (|Vila-Costas fc Edmunds] ll99ot ). 
The plot has a considerable scatter (~ 0.11 dex) around a 
mean correspo nding to the solar value, log(N/0)o = —0.86 
|Lodders!l2003h . and seems more compatible with a primary- 
plus-secondary origin for this element. The dispersion found 
in this relation is often explained as due to a time delay 
between the release of nitrogen and the release of oxygen 
back int o the interstellar medium by the respective stella r 
sources (|Vila-Costas fc Edmund3l993l : lMallerv et alj|2007h . 
Also H n regions in the galaxy M3 3 display an approxim ately 
constant value of the N/O ratio (|Bresolin et al.ll20101 ). with 
a mean log(N/0)~ —1-2, i.e. between the solar value and 
typical values for dwarf galaxies. The flat radial distribution 
of the N/O ratio is shown in the top panel of Fig. [7] 

The S/O ratio is constant as a function of 12+log(0/H), 
with a mean value log(S/0) = — 1.55 ± 0.07, consis- 
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Figure 5. Logarithm of (S+ + S ++ )/(0+ + 0++) versus + /0 
for our H II region sample (filled circles). The lines show the sulfur 
ionization correction factor (Eq. [4j for different values of a and 
assuming a constant value for log(S/0). 



the ICF given by the formu la (|Peimbert fc Costerd 1 19691 ; 
IStasinskalll975 lFrencrJll98lh : 



ICF(S+ + S ++ ) 



S+ + S H 



1-1- 



91 
o 



-l/a 



(4) 



where we have employed a = 2.5 (|Bresolin et al.l f2004h . 
Fig. 03 shows log[(S + +S ++ )/(0++0 ++ )] for our Hn re- 
gion sample versus the oxygen fractional ionization + /0. 
Eq. 2] has been overplotted for a — 2 and a — 2.5 and as- 
suming log(S/0) = — 1.5 and log(S/0) = — 1.6 (representa- 
tive values for our sample). Apart from three low-excitation 
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Figure 6. Trends in the abundance ratios relative to oxygen 
with 12+log(0/H). The horizont al dotted l ines i ndicate the corre- 
sponding solar value according to Loddcrs (2003). The filled circle 
co rresponds to the Hll region K932 with abundances measured 
bv lEsteban et al.l i2009l '). 
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Figure 7. Trends in the abundance ratios relative to oxygen 
with galactocentric radius. The horizontal dotted l i nes in dicate 
the corresponding solar value according to lLoddersi (2003). The 
filled circle c orresponds to the H I I region K932 with abundances 
measured bv lEsteban et al.1 (120091 ). 



tent, wi t hin e rrors, with the solar value (—1.50 ± 0.06, 
lLoddersi [2003). The Ar/O ratio is also essentially con- 
stant with oxygen abundance, but with a mean value of 
log(Ar/0)=-2.22 ± 0.06, slightly subsolar but marginally 



log^Ar/uj=— l.zz ± U.Ub, sngntly suosoiar out marginally 
compatible with the solar value C-2.14±0.09, lLodderdl2003M 



within the uncertainties. The Ne/O ratio, with a mean value 
of log(Ne/0) = —0.82 and a standard deviation of ~ 0.13 dex 
also shows a rather flat distribution in the whole range of 
oxygen abundance measured in our sampl e of Hll region s, 
and matches the solar value (—0.82 ± 0.11, Loddcrs 2003). 

Our data are compatible with a fiat distribution of the 
S/O, Ar/O and Ne/O ratios with galactocentric distance 
(see Figure [7|, but since our Hn region data are concen- 
trated at only two different galactocentric distances further 
data would be required to confirm these trends. 



4.4 He abundance 

Our Hn region spectra contain a number of He I recom- 
bination emission lines that can be used to estimate the 
He + abundance, as the emission line is proportional to the 
line emissivity and to the corresponding ionic abundance. 
Hen A4686 is not detected in any Hll region of the sam- 
ple and therefore the contribution of He ++ to the total He 
abundance can be neglected. 

The ionic helium abundance relative to hydrogen has 



been com puted from the predicted emissiv ities for the 
He I lines JPorter. Ferland. fc MacAdaml [20071 ) and for R/3 
|Beniamin et al.ll 1999h at the corresponding H n region elec- 
tron temperature, as computed in Sect. 14. ll 

The He 1 emission lines are affected by underlying 
stellar absorption, which needs to be co rrected for. We 
have followed the approach described by iKennicutt et al.l 
(2003). The corrections are based both on direct measure- 
ments of the He 1 EWs for O-type stars and He I EWs 
predicted by starburst models as a func tion of age. We 
have a pplied here the same corrections as IKennicutt et al.l 
(j2003T ). except for He 1 A4026 (0.6 A) and Hei A4922 
(0.3 A), which have been estimate d from more recent mod- 
els |Gonzalez-Delgado et al.ll2005T ). These corrections range 
from ~1% to 6% for the brightest lines (A5876 and A6678), 
-6-30% for A4922, and ~2-40% for A4471. 

The He + abundances computed from the brightest ob- 
served He 1 lines (A4026, A4471, A4922, A5876, A6678) are 
shown in Table [TOl together with the weighted average (col- 
umn 7) of all the individual values . The He + abundance s 
have also been computed using the iBeni amm et al.1 (|l999l ) 
emissivities. Differences in the a bundances derived using the 
more recent iPorter et all (|2007T ) emissivities are on average 
~2%, and in all cases below 5% and therefore below the 
He + /H + abundance uncertainties. 

The total helium abundances have been computed with 
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Table 10. Hc+ abundance. 



He+/H+ 
A4922 A5876 



ID 



A4026 



A4471 



A6678 



Average 0- 



2 

3 

17 

19 

25 

26 

33 

34 

35 



0.056 ±0.005 

0.086 ±0.010 

0.082 ±0.007 
0.083 ±0.010 



0.094 ±0.008 



0.063 ±0.006 
0.074 ±0.006 
0.081 ±0.008 
0.060 ±0.004 
0.083 ±0.007 
0.088 ±0.011 
0.058 ±0.006 
0.093 ±0.007 
0.087 ±0.008 



0.066 ±0.006 
0.076 ±0.006 
0.086 ±0.008 

0.085 ±0.007 
0.092 ±0.010 



0.069 ±0.010 
0.075 ±0.007 
0.094 ±0.011 
0.067 ±0.006 
0.088 ±0.008 
0.088 ±0.011 
0.054 ±0.006 
0.092 ±0.007 
0.084 ±0.010 



0.071 ±0.009 
0.079 ±0.007 
0.089 ±0.011 
0.068 ±0.006 
0.087 ±0.008 
0.088 ±0.010 
0.058 ±0.006 
0.081 ±0.007 
0.087 ±0.010 



0.066 ±0.004 
0.076 ±0.003 
0.088 ±0.005 
0.065 ±0.003 
0.086 ±0.004 
0.088 ±0.005 
0.057 ±0.003 
0.089 ± 0.004 
0.087 ±0.005 



a Weighted average of the He+/H + abundance ratio computed from all He I lines in the table 
(A4026, A4471, A4922, A5876, A6678). 
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Figure 8. Total He abundance versus galactoccntric distance 
(black filled circles) for our Hn region sample. The ionic He+/H+ 
abundances are plotted with open circles. The tot al and ionic He 
abund ances of the H II region K932 analysed by lEsteban et al.l 
are marked with a filled and open square, respectively. 



the analytical relation between the ICF( He) and the soft- 
ness parameter, T7=( 0+/0 2 +)/(S + / S 2+ ') dVflchez fc Page! 
Il988l), obtained by iBresolin et aT] ll2009al) from models 



of IStasinska. Schaerer. fc Leithererl l|200ll ). The ICF(He) 



ranges from 1.0 to 1.48 for our Hn region sample. The re- 
sulting total He abundances are shown in column 7 of Ta- 
bleland plotted versus gala ctocentric radius in F ig. El The 
figure includes region K932 |Esteban et al.|[2009t ). which is 
in very good agreement with our data for similar galacto- 
centric radius. The He abundance appears to be constant 
within our range of galactocentric distances, with an aver- 
age value He/H + = 0.092 ± 0.012, where the error represents 
the standard deviation between all regions. 



5 OXYGEN RADIAL ABUNDANCE 
GRADIENT 

The Hn region oxygen abundances in M31 obtained in this 
paper from the direct T e -based method are plotted as a 
function of galactocentric radius in Fig. [9] (blue squares). 
Th e linear regress i on to this data plus the H 11 region K932 
by lEsteban et al. (2009) yields the following solution: 

12+log(0/H) = 8.72(±0.18) - 0.028(±0.014) x R (kpc) (5) 

or equivalent ly, 



12+log(0/H) = 8.72(±0.18) - 0.56(±0.28) x (R/R25) (6) 

The slope that we obtain, —0.028 ±0.014 dex kpc -1 , repre- 
sents a rather shallow oxygen abundance gradient in M31. 
The central abund ance is only slight ly above solar (12 
±log(0/H) Q =8.69, lAsplund e t al. 2009). We are aware that 
our data are virtually representative of only two distinct 
galactocentric distances, and further data at different dis- 
tances from the M31 center would be desirable to better de- 
termine the abundance gradient from the T e -based method 
alone. 

This is the first time that an oxygen abundance gradient 
from auroral lines is obtained for M31. Previous attempts 
to measure the oxygen abundance gradient in this galaxy 
from H 11 regions were performed from diagnostics based on 
strong emission lines, because of the lack of the faint auroral 
lines necessary to measure the elect r on te mper atures. The 
early papers by Dcnncfe ld fc Kunthl ||l98ll ) and lBlair et al.l 
(1982) contain H 11 region samples of similar sizes (8 and 
11 regions, respectively) and make use of empirical ca libra- 
tions of the [O ill] /[Nil] line ratio (|Alloin et al.lll979l) and 
the R 23 parameter (R 23 =([0lll] + [Oll])/H/3, IPagel et~aT] 
Il979tl to derive the T e from which the O/H is estimated. 
Their O/H vs. galactocentric radius relations have a lot of 
scatter, but are in agreement with each other and seem com- 
patible with super-solar oxygen abundance in the M31 cen- 
ter (12±log(0/H) > 9.1) and an abundance gradient with 
a slope of approximately —0 . 03 de x kp c (from Fig . 4 an d 
Fig. 6. of lDennefeld fc Kunthl |l98ll ) and lBlair et all l|l982t l. 
respectively). This value is in agreement with our slope, but 
their metallicities are about 0.5 dex systematically higher 
than our determinations. 

Additional authors have recalculated the oxygen 
abundance gradient in M31, by applying different cal- 
ibr ations of the R23 parame ter to the data o btaine d 
by iDennefeld fc Kunthf (119811 ) and iBlair et all (|l982l ) 



Vila-Costas & Edmunds 



used the R23 calibration by 
Edmunds fc Pagell (Il984l). w i th adjustments at high and low 
metallicities |EdmundJl989l ; ISkillmanlll989l ) . They obtainec 
a somewhat steeper slope —0.043 dex kpc -1 . In contrast, 
IZaritskv et~aH l| 19941 ). using the same dataset, obtained a 
shallower abundance gradient (—0.018 ± 0.006 dex kpc -1 ) 
with a com bination of three different calibrations of the R23 
parameter llEdmunds fc Pagell 1 19841 ; iDopita fc Evans! Il986l ; 
iMcCall. Rvbski. fc Shields! 1 19851 ). In both cases a super- 
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Figure 9. Radial mctallic ity gradient fr om the direct T e -based method obtained from our work (filled blue squares). A direct abundance 
measurement from lEsteban et alj (l2009h for Hn region K932 is indicated w ith a cya n circle. A-, B- an d F-s upergiant oxygen a bundances 
from different authors are indicated with filled yellow l lTrundle et al.ll2002T) , green llVenn et al.ll2000h , re d iSmartt et al jfeoplh and pink 
llPrzvbilla et alj|2006l) stars. We also inlude the recent oxygen abundance determination for M31 PNe by lKwitter^T^lTil2012l ) with open 
dots. The solid straight line shows the best linear fit to the Ho region direct oxygen abundance versus galactocentric radius (see Sect.0, 
while t he dotted, dash ed a nd dot— dashed s t raight lines show the best linear fits to the oxygen abundances obtained from the IPilvuginI 
(200f ), Bresolin (2007) and Pcttini & Paeel (2004) strong-line methods, respectively, applied to the whole sample of Hll regions available 
in the literature (see Sect. UTTl l. 



solar central abundances (again 12+log(0/H) > 9.1) was 
obtained. 

The most recent spectroscopic observati ons of a sam- 
ple of Hll regions in M31 were performed bv lGalarza et all 
l| 19991 ). Their sample comprises 46 H II regions which the 
authors classify according to their morphological type. 
These authors inferred a metallicity gradient of —0.06 ± 
0.03 dex kpc" 1 from the dependence of the R23 parameter of 
the regions classified as center-brightened (~20 Hll regions) 
with the radial distance from the galaxy center. This gradi- 
ent is steeper compared to our and previous determinations 
of the oxygen abundance gradient of M31, but given its large 
uncertainty this result is still compatible with the shallower 
slope we find. Further spectrosco pic observations of si x H 11 
regions in M31 were obtained bv lBresolin et afl (|l999l ). but 
for a different purpose, and no abundance estimations were 
performed. 

More recent published values of the M31 oxygen abun- 
dance gradient from H 11 regions come from the reanalysis 
of the datasets described above dDennefeld fc Kunthlll98ll ; 
Blair et al ] | 19821 ; iGalarza et al.lll999T) using different empiri- 
cal calibrations dPaeel. Edmund~ fc Smith 1980: iMcGaughl 
ll99ll ; IZaritskv et allll994lPilvuginll200ll ). Depending on the 
empirical calibration applied, the slope ranges from —0.013 
to —0.027 dex kpc -1 , an d the central abundance varies be- 
tween 8.70 and 9.20 dex (|Smartt et all 1200 ll ; [Trundle et all 
120021 ). 



5.1 Strong-line methods 

The previous section makes it clear that the usage of dif- 
ferent calibrations of ratios of bright nebular emission lines 
(strong-line methods), which are readily accessible observa- 
tionally to determine the oxygen abundance gradient from 
a given data set, yields systematic and non-negligible dif- 
ferences on both the metallicity gradient and on the central 
(R=0 kpc) oxygen abundance. This is well known from the 
detai led analysis of different galaxies from many authors 
(e.g. iPerez-Montero fc Dfazl 120051; iKewlev fc Ellison! 120081 ; 
iBresolin et al .1 [2 009 al ibi: [Lopez-Sanchez fc Estebanll2010l ). In 
particular, strong- line methods calibrated via grids of pho- 
toionization models generally provide higher metallicities 
than those calibrated empirically via sets of H 11 region oxy- 
gen abundance determinations obtained from T e estimates 
from detected auroral lines. 

Given the still limited number of regions for which we 
have a direct determination of the electron temperature in 
M31, it is useful to apply different bright-line metallicity 
indicators to the whole set of M31 H 11 regions observed 
so far to enlarge the range in galactocentric distance and 
the number of regions for a better statistical significance 
of the gradient parameters. As a first step, we have ap- 
plied strong-line methods to the whole set of 31 Hll re- 
gions analyzed in this paper (those with and without de- 
tections of auroral lines), limiting our choice to six widely 
used m ethods: (a) The R23 met hod, as calibrated empiri- 
cally by IPilvuginI (|200ll . POl) and lPilvugin fc Thuanl (|2005l . 




Figure 10. Radial oxygen abund ance grad i ent i n M31 from differe nt calibrations of st rong-line ratios applied to the Hn regions analyzed 
in this paper (filled symbols): (a) iPilyugin (20 0l], P01) and IPilyugin &; Thuanl J2005I . PT05) empirical calibr ations of the R23 indicator. 
(b) R 23 strong-line me thod as calibrated by IMcGaughl <1991f) . (c) N202 (=[N III /[O III) as calibrated by Bresolinl ( |2007| . B07) and 
IKewlev fe Dopital ||2002| . KD02) (d) N2 =log([N Il]6583/Ha) as empirically calibrated by iPettini fc Paeej ^2004 PP041. The dotted 
straight lines show the best linear fit to the oxygen abundance gradient obtained from the corresponding strong-line methods. The fit 
parameters are shown in Table [TT] The solid straight line shows the best fit to the oxygen abundances obtained from direct measurements 
of T e (Eq. |5}. Filled squares indicate H II regions with detections of auroral lines, while filled hexagons indicate regions with no detections 
of auroral lines. The open symbols show the direct T e — based oxygen abundances derived in this paper. 



PT05) from spectroscopic data of H 11 regions with measured 
T e -based oxygen abundances, (b) The R2 3 method, ascal- 
ibrated from theore tical model grids by iMcGaughl fl99ll. 
M91), through the iKuzio de Narav. McGaugh, fc de Blokl 
l|2004 ) parametrization, which takes into account the ioniza- 
tion parameter of the H n regions, (c) The N202 method, 
i.e. the empir ical calib r ation of the [N 11] A6583 /[On] A3727 
line ratio by iBresolinl (|2007l . B07) with extragalactic Hn 
regions with T* e -bas ed oxygen abundances, and the cali- 
bration obtained by IKewlev fc Dopital (|2002l ) from stellar 
population synthesis and photoionization models, (d) The 
N2= log([N h1A 6583/Hq) method, empirically calibrated by 
IPettini fc Pagell (|2004 PP04) with Hn regions with Te- 
based oxygen abundances or detailed photoionization mod- 
elling. 

The results from the application of these methods to our 
Hll region sample are plotted in Fig. HOI as a function of the 
galactocentric radius. For the three R,23-based calibrations 
we have used the upper branch, as indicated by the discrim- 
inator [Nii1A6583/[Oii 1A3727 (|Kuzio de Narav et all 12004 
IKewlev fc Ellison] 2008). For comparison we also plot (with 
open symbols) the Te-based oxygen abundances derived in 
this paper together with the best linear fit to the data (solid 
line). The oxygen abundances obtained with the R23 param- 



eter as calibrated bv lMcGaughl (Il99ll) and the ones obtained 
with the N202 calibration of IKewlev fc Dopital (|2002l ) are 
~0.4-0.7 dex larger than the ones obtained from the other 
methods, as expected when comparing empirical and the- 
oretical model grid calibrations. Table [TT] shows the best 
linear fit parameters for the oxygen gradient obtained from 
the different strong-line methods. 

It is apparent from Fig.[l0]and Table [TTI that the scatter 
varies between the different O/H indicators, being larger for 
the R23 -based methods and N2, and smaller in the case 
of N202 which presents the highest correlation coefficient 
(around —0.85 for both calibrations of N202, with a root 
mean squared (rms) residual scatter from the fit of 0.05 dex 
and 0.03 dex for the IBresolinl l|2007l ) and IKewlev fc Dopital 
(2002) calibrations, respectively (cf. 0.12 dex for N2 and for 
the PT05 calibr ations and 0.14 dex for R23 as calibrated by 
lMcGaugblfl99ll) . 

Fig. [10] and Table [TTJ show that the s lope of the oxy - 
gen abundance gradient obtained from the IPilyugin! (l200ll ) 
calibration of the R23 par a meter , from N2, and from N202 
as calibrated by IBresolinl |2007T ) agree within the errors, 
and also agree with the oxygen abundance gradient ob- 
tained from the direct T e -based ox y gen a bundances ob- 
tained in this paper. The IMcGaughl (|l99lh calibration of 
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the R23 parameter agrees marginally with the N2 calibra- 
tion and with the T e -b ased oxygen grad i ent. T he N202 
method as calibrated by IKewlev fc Dopital (|2002l ) and the 
IPilvugin fc ThuanI (2005) calibration give shallower oxygen 
abundance gradients , but we note that the cor relation of the 
abundances from the lPilvugin fc Th uan ( 2005]) method with 
galactocentric radius is poor (correlation coefficient ~ —0.2). 

The six strong-line methods have also been applied to 
all spectroscop ic data available in the lit e rature for Hu re- 
gions in M31 llDennefeld fc Kunthl Il98ll; iBlair et al.lll982l; 



Galarza et al. I 1 19991 ; iBresolin et al. I 1 19991 ; Estc ban et al.l 
20091 ). The galactocentric distances have been recalculated 



from their corresponding equatorial coordinates (as reported 
by the various authors) and using the same galactic pa- 
rameters as for the estimation of the galactocentric dis- 
tances of our Hn region sample (see Table |5J). Only the 
Hu regio ns morphologically cl assified as center-brightened 
from the iGalarza et al.l l|l999h sample hav e been consid- 
ered. We also note that these authors and IBresolin et al.l 
i|l999l ) only provide the total [Nil] (A6548+A6583) emission. 
In order to apply the N2 and the N202 methods to these 
samples, we have estimated the expected [Nll]A6583 emis- 
sion by assuming [Nn]A6583 = 3x [N ll]A6548. Also, in the 
cases in which [O lll]A4959 is not detected, we assumed that 
[Om]AA4959+5007=(4/3) [Om]A5007. The full Hu region 
sample comprises 85 objects. 

The results are plotted in Fig. [TT] and the correspond- 
ing least-square fits parameter shown in Table 1111 We ob- 
tain a similar behaviour as for our data alone: we observe 
an agreement between the T e -based oxygen gradient and the 
slopes obtained from strong-line methods with the exception 
of the IPilvugin fc ThuanI ( 20051) empirical c a librat ion and 
the N202 as calibrated by Kewlev fc Dopital l|2002l ). These 
methods yield a substantially shallower slope (— 0.009±0.005 
an d —0.011 ± 0.002 dex kpc ~ , respectively). In the case of 
the lPilvugin fc ThuanI (2005) calibration, we find again that 
the correlation is poor ( correlation coefficient —0. 24). This 
lack of agreement of the lPilvugin fc ThuanI (|2005l ) calibra- 
tion with the r e -based o xygen abundances is surpris i ng an d 
it is in contradiction with Lope z-Sanchez & Estcban (2010), 
who conclude that the lPilvugin fc Thuan 1 20051 ) calibration 
is nowadays the best empirical method to estimate oxygen 
abundances when auroral lines are not detected. 

We should note that in all empirical calibrations we 
slightly i mprove the correlatio n coefficient and the scatter 
when the lGalarza et all (|l999l ) sample is not taken into ac- 
count, but the slope and intercept of the fits remain un- 
changed. 

We have computed the weighted average of the slope of 
the oxygen abundance gradient from the strong-line meth- 
ods with lower scatter and with a better match with our 
direct T e -based oxygen abundances (i .e. with the N202 
m ethod as calibrated by IBre solin 2007 ], the N 2 calibrated 
bv lPettini fc Pagellliool and the lPilvuginll200ll calibration) 
with weights equal to the reciprocal of the errors, yielding 
the value 



dlog{0/H)/dR = -0.023 ± 0.002 dex kpc" 1 



(7) 



that we adopt as the slope of the oxygen abundance gradient 
in M31. 



6 COMPARISON WITH OTHER 
METALLICITY INDICATORS 

We have compared the H n region abundance gradient with 
the results obtained from young, BA-type supergiants and 
from planetary nebulae, which trace an older stellar pop- 
ulation. Despite the relatively small distance to the An- 
dromeda galaxy and the importance of this extragalactic 
neighbour in the context of galactic evolution studies, the 
analysis of a sizeable sample of blue supergiant spectra 
for the determination of stellar metallicities is still lack- 
ing. The situation has remained virtually unchanged dur- 
ing the past decade, since the public ation of stellar metal- 
licities for 7 B-type supergiants by iTrundle et aTJ l|2002l ). 
and, somewhat paradoxically, lags behind the quantitative 
spectroscopic study o f massive stars in considera bly more 
distant galaxies (e.g. iKudritzki et al.ll2008l . |2012|). One of 
these stars had been originally analyzed bv ISmartt et al.l 
(2001). Additionally, t wo A- and o ne F -type supergiants 
were investigated by Venn et alj (|2000| ). More recently 
IPrzvbilla. Butler, fc Kiidritzkil ( 20061 ) revisited the metallic- 
ity derived for one of the A-type stars with more accurate 
non-LTE calculations, and analyzed a new A-type star. It 
sh ould be pointed out t hat the B-type metallicities obtained 
bv lTrundle et all (2002) were based on an LTE analysis, ex- 
cept for one case, in which the non-LTE metallicities were 
found to be consistent, within the relatively large (between 
~0.2 and 0.4 dex) errors, even though they are 0.07 dex lower 
in the case of oxygen (a similar result is obtained comparing 
the LTE vs. non-L TE analysis of the star in common with 
ISmartt et aT]|200ll ). 

These studies highlighted the difficulty of comparing the 
stellar and the nebular metallicities, due to the systematic 
differences in Hn region abundances obtained from differ- 
ent strong-line methods. For this reason, the results of the 
comparison were inconclusive. Our new study of H n regions 
in M31 was mostly motivated by the need to acquire direct 
Te-based nebular abundance. As this paper shows, we have 
partially succeeded in our goal, but the detection of auroral 
lines in H n region spectra remains elusive for M31 even with 
lOm-class telescopes. 

In Fig.[9]we include the O/H abundances derived for the 
B-type supergiant stars by Trundle et al. (2002, yellow star 
symbols) and Smartt et al. (2001, red star symbol), and for 
the AF stars studied by Venn et al. (2000, green star sym- 
bols) and Przybilla et al. (2006, pink star symbols). When 
available, only the non-LTE result is shown. Galactocentric 
distances have been recalculated from our adopted M31 pa- 
rameters (see Table |2|). The comparison is somewhat com- 
plicated by the large errors in the B supergiant abundances, 
but it is obvious that a systematic offset is present in relation 
to the direct abundances we determined from Hn regions. If 
we focus on the A-type supergiants (galactocentric distance 
~11.7 kpc), which have more accurate abundances, the neb- 
ular O/H ratio lies ~0.3 dex below the stellar value. On the 
other hand, among the gradients obtained from strong-line 
diagnostics, we find t hat the N2 method, as calibrated by 
IPettini fc Pagell J 20041 ) . is in acceptable agreement with the 
absolute value provided by the stars (the difference relative 
to the stars is reduced to ~0.1 dex). This situation contrasts 
the results obtained by our group in other, lower metallic- 
ity galaxies, for example NGC 300, where excellent agree- 
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R (kpc) R (kpc) 

Figu re 11. Same as Fig. 1101 but for the H II regions ana l yzed i n this paper plus H II regions spectroscopically observed by other authors: 
IDennefeld fc Kunthl jl98lll. yellow triang les; iBlair et al.1 1 19821) . red diamonds; H II regions classified as ' center-bright' bvlGalarza et all 
(1999), green stars; iBresolin et al.l lll999l ). pink inverted triangles; region K932 (cyan circle) observed by Estcban ct al. (2009j), and Hn 
regions analyzed in this paper (blue squares and hexagons, as in Fig. HOP . Open symbols indicate Hn regions with uncertain values in 
the relevant line ratios for a given empirical calibration, as stated by the corresponding authors. The dotted straight lines show the best 
linear fit to the oxygen abundance gradient obtained from the strong-line methods applied to the whole Hn region dataset (excluding 
regions with uncertain line ratios). The solid straight line shows the oxygen abundance gradient obtained from the T e — based method 
(Eq.0. 



Table 11. Linear least-squares fit parameters to the oxygen abundance distribution across 
the M31 disk obtained from different strong-line methods (see Sect. 15. Il l applied to our 
Hn region sample and to the whole set o f spec t roscopic data on M31 H II regions pub - 
lished so far: this pap e r. IDennefeld & Kunthl lll98ll) . lBlair et al.l l ll982l) , lGalarza et al.l fl999l) . 
IBresolin etafl \ 19991) . lEsteban et al.l J200S| ). 



Method 


slope 


zero point 


Correlation 


No. 


Scatter^ 




(dex kpc ) 


(dcx) 


coefficient 


data points 


(dcx)" 


Our H II region sample 


R.23 (P01) 


-0.021 ±0.004 


8.82 ±0.05 


-0.69 


25 


0.08 


R23 (PT05) 


-0.007 ±0.007 


8.45 ± 0.08 


-0.20 


25 


0.12 


R23 (M91) 


-0.037 ±0.008 


9.39 ± 0.09 


-0.70 


25 


0.14 


N202 (B07) 


-0.024 ±0.003 


8.77 ±0.03 


-0.86 


28 


0.05 


N202 (KD02) 


-0.015 ±0.002 


9.43 ± 0.02 


-0.84 


28 


0.03 


N2 (PP04) 


-0.027 ±0.006 


8.88 ±0.07 


-0.65 


31 


0.12 




H11 


regions from all authors 






R23 (P01) 


-0.020 ±0.004 


8.76 ±0.05 


-0.55 


70 


0.12 


R23 (PT05) 


-0.009 ±0.005 


8.45 ± 0.06 


-0.24 


70 


0.15 


R23 (M91) 


-0.029 ±0.005 


9.25 ±0.06 


-0.59 


70 


0.15 


N202 (B07) 


-0.023 ±0.002 


8.74 ± 0.03 


-0.75 


75 


0.08 


N202 (KD02) 


-0.011 ±0.002 


9.39 ±0.02 


-0.66 


75 


0.05 


N2 (PP04) 


-0.024 ±0.003 


8.87 ±0.04 


-0.64 


82 


0.11 



a rms residual scatter, estimated as the square root of the arithmetic mean value of the 
squares of the difference between fitted and measured data. 
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ment is found between the metallicities obtained for stars 
and Hn regions using both the direct T^-based method and 
N2 (|Bresolin et al.ll2009ah . A situation si milar to what we 
encounter in M31 is also found for M81 (|Patterson et al.l 
l2012l ; lKudritzki et alj20l3 ) . A possible interpretation of this 
discrepancies will be provided in t he next section. 

Recentlv lKwitter et al.l ()2012h published the abundance 
analysis of 16 PNe from [Oin]A4363 detections in the outer 
disk (R > 20 kpc) of M31. The radial distribution of the 
O/H ratio is included in Fig. [5] (open circles; we recalcu- 
lated the galactocentric distances from our adopted orien- 
tation parameters, Tabled finding only minor differences). 
We note a striking offset between Hn regions and PNe di- 
rect abundances (although we note that their target 16 falls 
near the Hn region locus). It is not clear whether the two 
samples of objects are directly comparable, however, due to 
the fact that the PNe belong to the outer disk of the galaxy. 
Difficulties in the interpretation of the offset might also arise 
from the complicated chemical history of M31 (e.g. Davidge 
ct al. 2012 for a recent summary). Finally, w e note that the 
slope of the PNe O/H gradient obtained bv iKwitter et all 
(|2012l l. -0.011 ± 0.004 dexkpc" 1 , is somewhat flatter than 
the one we derive from Hn regions. We add that a flatter 
gradient for PNe compared to H n regions has also been re- 
cently determined in the galaxy NGC 300 (Stasinska et al. 
2012, in preparation), and corrobora tes recent model predic 
tions of galactic chemical evolution ijPilkington et al 



predic- 
I l2012h . 



7 DISCUSSION AND CONCLUSIONS 

One of the main conclusions obtained in this paper is that 
the slope of the oxygen abundance gradient in M31 from 
the analysis of the Hn regions is quite robust, approxi- 
mately —0.023 dexkpc^ 1 , regardless of the choice of abun- 
dance diagnostic used (wit h the notable exception of the 
IPilvurin fc Thuanll2005l and lKewlev fc Dopitall2002l calibra- 
tions). However, the ~0.3 dex offset between the direct abun- 
dances (obtained from the auroral line detections) and the 
stellar abundances is somewhat disturbing, especially con- 
sidering the fact that the direct abundances also lie sig- 
nificantly below those obtained from strong-line methods 
that, by construction, should be in agreement with the di- 
rect determinations. Since the discrepancy could be a result 
of the higher metallicity of M31 compared to other galaxies 
where such offset is not observed, we suggest an interpreta- 
tion based on biases affecting the direct abundances. As is 
well known, with increasing O /H the detection of the auroral 
lines in the spectra of H n regions becomes increasingly dif- 
ficult. This is due to the exponential dependence of the line 
emissivities on the electron temperature, and the decreas- 
ing nebular temperatures with increasing metallicity due to 
more efficient cooling. As a result, very few auroral line de- 
tections exist in the literature for extragalactic Hn regions 
with metallicity above solar. Therefore, given the metallic- 
ity (around the solar value) found for the supergiant stars in 
M31 it is not surprising that we only have a few [O ill] A4363 
detections. Under these conditions we would preferentially 
detect this line just for some of the hottest (less metal en- 
riched) H n regions, which leads to a bias in the derived av- 
erage oxygen abundance: the abundance estimated from the 
few [Om]A4363 detections would underestimate the mean 



O/H ratio. This would explain why the auroral line-based 
abundances lie below the strong line abundances obtained 
from methods such as N2 and N202, which have been cali- 
brated using H II regions with detections of auroral lines. If 
we assume an observational abundance scatter of ~ 0.07 dex 
(comparable to what has been measured in other nearby 
galaxies, see 



Brcsolin 



2011 



;o wnat nas p een mcasu rea m otner n earby 
Bresolin et all l2009al . iKennicutt et al l 12003 . 



lj), we can expect to systematically underes- 
timate the oxygen abundance by as much as ~0.15 dex 
if we assume that our auroral line detections refer to the 
hottest nebulae only, i.e. are not representative of the mean 
Hn region population. Most of the remaining discrepancy 
with the stellar abundances could be explained by depletion 
of oxygen onto dust grains, which ha s been estimated to 
be ~0-12 dex at high metallicity (e.g. iMesa-Delgado et al.l 
120091 . iPeimbert fc Peimbertl [2010 1 . If we account for these 
effects, the central abundance of M31 would be approxi- 
mately 12+log(0/H) = 8.94 ± 0.10 dex; i.e. between 1.4 and 
2.2 times the solar value. 

Finally, we also point out that a discrepancy between 
direct nebular abundances and stellar abundances increasing 
with metallicity would naturally result from the presence of 
inhomogen eities in the tempera ture distribution within an 
Hn region (Pcimbcrt ct al. 20 0^) or by the effect recently re- 
ported bv lNicholls. Dopita. fc Sutherland! j2012), who pro- 
posed that the energy distribution of electrons in H II regions 
deviates from a Maxwell-Boltzmann one. Further tests of the 
mechanisms proposed here to explain the systematic offset 
between T e -based Hn region abundances and stellar metal- 
licities would require to increase the sample of Hn regions 
with direct abundance determinations for metal-rich (super- 
solar) galaxies, together with a program aimed at measuring 
metallicities of blue supergiant stars. In the specific case of 
M31, in particular, more accurate stellar abundance data for 
a large sample of B-type stars would be highly desirable. 
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